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Abstract 

In this paper we describe how an idea centered on the concept of self-saturation 
allows several improvements in the computation of Gröbner bases via Buchberger's 
Algorithm. In a nutshell, the idea is to extend the advantages of computing with 
homogeneous polynomials or vectors to the general case. When the input data 
are not homogeneous, we use a technique described in Section 2: the main tooi 
is the procedure of a self-saturating Buchberger's Algorithm, and the main result 
is described in Theorem 21. Another strictly related topic is treated in Section 3 
where a mathematical foundation is given to the sugar trick which is nowadays 
widely used in most of the implementations of Buchberger's Algorithm. A special 
emphasis is given in Section 4 to the case of a single grading, and Section 5 exhibits 
some timings and indicators showing the practical merits of our approach. 
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Introduction 



Starting from the sixties, when implementations of Buchberger's famous al- 
gorithm for computing Gröbner bases became practically feasible, it has been 
clear that mainly three of its steps can be optimized. They are the minimal- 
ization of the set of critical pairs (see for instance (0)), the optimization of 
the reduction procedure (see for instance (0)), and the sorting used to process 
the critical pairs during the algorithm. The last aspect is less important when 
the input polynomials or vectors are homogeneous and the algorithm proceeds 
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with an increasing degree strategy. But what happens if the input data are 
not homogeneous? 

A first answer to this question was given in the late eighties. It prescribed to 
homogenize the input data, run the algorithm, and then dehomogenize the 
computed Gröbner basis. This strategy is indeed quite simple and in many 
cases works fine. lts big advantage is that critical pairs are sorted by increasing 
degree and after a degree is completed the algorithm never goes back to it. 
The disadvantage is that often it computes too large a set of polynomials or 
vectors. 



Quite soon (we are speaking of the early nineties) a new tooi entered the 
game, the sugar strategy (see (lid ) and Section[3]). In a nutshell, the idea was 
to keep the data non homogeneous, but process the critical pairs as if they 
were coming from true homogeneous data. This goal is achieved with the 
help of a manipulated degree called sugar which substitutes the true degree. 
Although a complete theoretical background was not laid out, the idea gained 
popularity. Not much later, paper (Il4l ) described an implementation in the 
computer algebra Bergman which uses a way to improve the ordering of the 
critical pairs. Also that source was lacking a solid theoretical foundation and 
it did not gain the same popularity as the sugar strategy. 



Recently, inspired by the new development of CbCbA which will lead to the 
long awaited CbCbA 5, we decided to explore some features of Buchberger's 
Algorithm in great detail. The main purpose was to give a solid theoretical 
background to both the sugar strategy and the strategy of selection of critical 
pairs. We believe that we achieved both goals, so let us explain how. After 
recalling more or less well-known facts about the homogenization process in 
Section [T], we move quickly to the construction of what is called the Self- 
Saturating Buchberger's Algorithm. 



To do that, in Section [2] we prove several properties of the saturation (see 
Proposition[9]), define new notions such as a-SatGBasis and a-DehomBasis , 
and prove Theorem [13] where all these notions are fully compared. With 
the aid of this result we define and study several variants of Buchberger's 
Algorithm, the Weak Self-Saturating Buchberger's Algorithm and the Self- 
Saturating Buchberger's Algorithm, and finally prove the desired main result, 
Theorem [T9l It simply says that the computation of a Gröbner basis when the 
input data are inhomogeneous, can be performed by running any Weak Self- 
Saturating Buchberger's Algorithm. The inspiration to achieve this goal came 
not only from the above mentioned paper (14), but also from the paper (jïl) 
where similar strategies were described for the efficiënt computation of toric 
ideals. 



It is also noteworthy to mention the fact that the variants of the Weak Self- 
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Saturating Buchberger's Algorithm include the usual Buchberger's Algorithm 
as well as the algorithm obtained by homogenizing the input data, run the 
algorithm, and then dehomogenizing the computed Gröbner basis. 

Section [3] is devoted to give a solid foundation to the sugar strategy which, as 
we said, is already used in several computer algebra systems. To describe it in 
joking mode we could say that the idea is to make a recipe by adding some 
sugar to the degree of the inhomogeneous vectors and make them sweeter in 
Buchberger's Algorithm. The main result is Proposition[25] which describes the 
behavior of the sugar during the execution of every variant of Buchberger's 
Algorithm introduced in the previous section. With this result we can combine 
the tools of Section [2] with the sugar strategy. 

Section H] treats the case of a single grading and shows how in that situation 
better results can be achieved (see Theorem [271 and its coroUaries) . The cur- 
rent implementation in CbCbA deals only with the case of the single gradings, 
shortly to be extended to the general case, and the final Section [5] shows its 
excellent behavior on a selected bunch of examples. 

Of course we are aware of many algorithms which optimize the computation 
of some Gröbner bases simply by going around the problem. Among many 
others we could recall the Gröbner walk algorithm, the FGLM algorithm. 
But we want to make it clear that our goal here is to optimize Buchberger's 
Algorithm, not to find alternative strategies to compute Gröbner bases. 

As a side remark we observe that every Self-Saturating Buchberger's Algo- 
rithm is fuUy compatible with the SlimGB strategies developed in (j2|) and with 



the Hilbert driven algorithms (see (Il3l ) and (p)). The integration and inter- 



play of these approaches will be the subject of future work. Finally, the readers 



should know that the basic terminology is taken from the two books (jlll ). (jl2l ). 



1 Preliminaries 



We assume the basic terminology and facts explained in (jl2l ). Section 4.3 and 
Tutorial 49. Some of them are explicitly recalled for the sake of completeness, 
hence most of the section contains either well-known facts or easy generaliza- 
tions of well-known facts. 



1.1 Homogenization in a polynomial ring 



In this subsection we generalize the natural concept of homogenization to the 
multigraded case. 
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Wc let K hc a field and P = K[xi, . . . ,Xn] a polynomial ring. Then we 
take a matrix W G Matm,n(^) of rank m > 1 and new indeterminates 
Ui^-'-iUm called homogenizing indeterminates. Moreover, we equip the 
polynomial ring P = K[yi, . . . , y^, xi, . . . , Xn] with the grading defined by 
the matrix W — {lm \ W), where lm denotes the identity matrix of size m . 

Given m-tuples of integers Vj = [aij, . . . ,amj) , j = we consider 

the tuple (ci, . . . , c^) where = max{afei, . . . , aks} , k = 1, . . . , m, and call 
it Top(t'i, ...,Vs). 

Definition 1 Let ƒ G P \ {0} and F e P. 

(1) Write ƒ = Ci^i + • • • + Cstg with ci, . . . , G \ {0} and distinct terms 
tl, ... , tg G T"^. Then the tuple Top(deg^(ti), . . . , deg^y(tm)) is called the 
top degree of ƒ with respect to the grading given by W and is denoted 
by TopDeg^(/). 

(2) For every j' = 1, . . . , s, we let deg^y(tj) = (ry, . . . , Tmj) G Z"^ and let 
(/xi, . . . , /i^) = TopDeg^(/). The homogenization of ƒ with respect 
to the grading given by W is the polynomial 

For the zero polynomial, we set 0*^°™ = 0. 

(3) The polynomial F'^^^ = P(l, . . . ,l,xi, . . . , Xn) G P is called the deho- 
mogenization of F with respect to yi,...,ym- 

Given an ordering r on T", the monoid of power- products in P, we want 
to extend it to T"^+", the monoid of power-products of the homogenization 
ring P . 

Definition 2 We consider a monoid ordering r on T", and the relation 
on T™"*"" which is defined by the foUowing rule. Given two terms ii, t2 £ T"^^"' , 
we say that ti >^w t2 if either 

deg^ih) > deg^(t2) 

or 

deg^(ti) = deg^(t2) and tf^ t^^^ 

We call the extension of r by ly. If it is clear which grading we are 
consider ing, we shall simply denote it by r. 

We recall that the grading represented by the matrix W is said to be positive 
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if each column of W has some non-zero entry and the first non- zero entry is 
positive. 

Proposition 3 Let t be a monoid ordering on T" and r its extension by W. 

(1) The relation T is a deg-^^ -compatible monoid ordering on T™"'""'. 

(2) If W is positive, the relation t is a term ordering on Y^^^ . 

(3) Let F & P be a non-zero homogeneous polynomial. Then there exist 
si, . . . , s„ G N such that LT^(F) = • ■ ■ y^"^ ■ LT^{F^^'') . 

(4) If T is of the form t = Ord{V) for a non-singular matrix V G Mat„(Z), 
then we have T = Ord (^^^ . 



Proof. For the easy proof see (jl2l ). Proposition 4.3.14 and Lemma 4.3.16. □ 



Remark 4 If r is deg^;[/ -compatible then r = Ord . Therefore we have 



_ lm W \ I lm. W 

T = Ord ( O w ] = Ord -im o 

O V' I \ O V' 



If m = 1 it follows that r is of ï/i-DegRev type (see (jl2l ). Section 4.4) with re- 
spect to deg^y . In particular, if m = 1 and r = DegRevLex where Deg denotes 
the Standard grading on P it is more common to write P = K[xi, . . . ,Xn,y] 
with the homogenizing indeterminate at the end, then we have r = DegRevLex 
where Deg denotes the Standard grading on P. 



1.2 Homogenization in a free P-module 



In this subsection we generalize the multihomogenization procedure to the 
case of free modules. 

Let r be a positive integer, let F denote the free P-module P^ and let 
Cl, . . . , be the vectors of the canonical basis of F . Then let 5i, . . . , 5^ G Z"^ 
and let F be the graded free P-module F = 0[=i P{—6i) where the de- 
grees of ci, . . . ,er are 6i, . . . ,6r respectively. We denote by T"(ei, . . . , Cr) the 
monomodule made by the terms t-Ci e F with t G T" , and by T"^~^"'{ei, . . . , e^) 
the monomodule made by the terms t ■ Ci E F with t G T™+" . Henceforth, 
when we consider module orderings on T"'(ei, . . . , e^) we always mean mod- 



ule orderings which are compatible with a monoid ordering on T" (see ( llll ) 
Definition 1.4.17). The following definition is a natural generalization of Def- 
inition [2l 



Definition 5 We consider a module ordering a on T"(ei, . . . , e^) and the 
relation on T™"'""(ei, . . . , e^) which is defined by the following rule. Given 
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tiCi, t2ej G T™'+"'(ei, . . . , e,.) , we say that tiCj >-^w t2ej if either 

deg^^(tiei) > deg^{t2ej) 

or 

degy^itiCi) = deg^{t2ej) and tf '"ci >a tg^'^ej 

We call the extension of o" by VT. If it is clear which grading we are 
considering, we shall simply denote it by W. 

Proposition 6 Let a be a module ordering on T"(ei, . . . , e^) , and let a be 

its extension by W. 

(1) The relationa is a deg-y^ -compatible module ordering on T™'+"(ei, . . . , e,-) . 

(2) If W is positive, then a is a module term ordering on T^~^^{ei, . . . , e^) . 

(3) Let U E F be a homogeneous non-zero vector. Then there exist non 
negative integers Si, . . . , such that UY-^{U) = yl^ ■ ■ -y^ ■ \jTa{U'^^^) . 

Proof. It is an easy generalization of Proposition [31 □ 

Analogously to Definition [T] one defines the homogenization and dehomoge- 
nization of vectors, and with the following proposition we recall some easy 
results about homogenization and dehomogenization we will need to prove 
Theorem [T3l 

Proposition 7 Let M be a submodule of F which is generated by vectors 
Vi, . . . ,Vs, and let N be a graded submodule of F which is generated by ho- 
mogeneous vectors Vi, . . . ,Vt. 

(1) We have [M'^^'^Y^'^ = M. 

(2) The homogenization of M can be computed via the formula 



M^°- = «-,..., yf-) :_ {yi---y^r 

(3) The dehomogenization of N can be computed via the formula 



Proof. It is an obvious generalization of (jl2l ). CoroUaries 4.3.5.a and 4.3.8. □ 



2 Self-Saturating Buchberger's Algorithm 



This section starts with some properties of the saturation and continues with 
the proof of the main facts (see Theorem [T3l) which will eventually lead to 
the algorithm for computing inhomogeneous Gröbner bases. After recalling 
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the definition of a remainder, we write the body of Buchberger's Algorithm to 
help the reader spotting the differences when we describe some of its variants 
(see Theorem [Tül) . The section ends with the main Theorem[2Il We keep the 
notation introduced before, in particular, we let cr be a module ordering on 
T"(ei, . . . , Cr) , and let a be its extension by W. 

2. 1 Saturation 



Definition 8 Let U G -F be a homogeneous vector. We denote (^[/^eh^hom 
by f/^^* and we call it the saturation of U . Let be a graded submodule 
of F. We denote (iV^-^h^hom jysat ^-^^ saturation of A^. 

We are going to illustrate some fundamental properties of the saturation. 
First, it is useful to recall Proposition Oc where we showed that there exist 
si, . . . , s„ e N such that the formula LT^{U) = ■■■yf^ ■ LT^(f/'i^^) holds 
true. 

Proposition 9 (Properties of the Saturation) Let vi,...,Vs be vectors 
in F , let U E F he a homogeneous non-zero vector. 

(1 ) There exist ri, . . . , G N such that LT^(f/) = y[' ■ ■ ■ y''^ ■ ISI^iW^^) . 

(2) Comparing (1) with the formula Ul-^iU) = yl^ ■■■y^ • l^aiU'^'^^), we 
have Ti < Si for i = 1, . . . , m . 

(3) We have (LT^(f/))'^^t = LT^{U'^''^). 

(4) If N is a graded submodule of F, then we have the following equality 
N^-t = N ■._{yi---ymr. 

(5) If N is a graded submodule of F , then we have the following equality 
jysat ^ ^^sat | y ^ jy^ y fiomogcncous) . 

(6) If N is a graded submodule of F , then we have the following equality 

sat ^ deh jy^eh 

(7) If M = {vi, . . . ,Vs) is a submodule of F, then we have the following 
equaUty M^"'" = . . . , y^""^)'^' . 

Proof. Condition (1) follows from the definition. To prove (2) we denote by (*) 
the formula hT^iU) = y{^ ■ ■ ■ |/^™ ■ LT„(f/'^^'^) in Proposition [1(3). We observe 
that {U^^^Y'^^ = Lf^^^, hence, if we apply (*) to W'^^ we get the equality 

LT-s^(U^^^) = yl^ ■ ■ -y^ ■ UT^{U^^^) for suitable natural numbers s[, . . . , s^. 
Using (1) and (*) we get + s- = Si for i = 1, . . . ,m. Condition (3) fol- 
lows from condition (2). Next we prove (4). Let Vi, . . . ,Vt be homogeneous 
vectors which generate N . Using Proposition [71 we deduce the following 
equality A^*^^* = {Vl^^ , . . . ,V^^^) :-p {yi - ■ ■ym)°^ ■ It remains to show that 
-F (^1 ■ ■ ■ l/m)°° = (Vf^*^*, . . . , V/*^*) ■ ■ ■ ym)°° • The inclusion C is a con- 
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sequence of the obvious relation N C {Vi'^^, . . . , V^*^*^*) , while the inclusion D 
foUows from the observation that V^^^ G :-p (ï/i ■ ■ ■ym)°° for i = 1, . . . ,t. 
Condition (5) foUows from the definition. Clearly (6) follows from (4) and 
finally, to prove (7) it sufïices to combine (4) with Proposition [71(2). □ 

Definition 10 Let be a graded submodule of F and let Vi, . . . ,Vt G 
be non-zero homogeneous vectors. 

(1) The set {Vi, . . . ,Vt} is called a W-SatGBasis for if it is a a-Gröbner 
basis of a graded submodule N of F such that N^^^ = N^^^ . 

(2) The set {Vi, . . . , 14} is called a a -DehomBasis for A^ if {Vf"'^, 1/^'^"^} is 
a cr-Gröbner basis of A^'^'^'^. 

Proposition 11 Let M = {vi, . . . ,Vs) be a suhmodule of F, denote by N 
the module {v\°'^, v^""^) , and let {Vi, . . . , 14} be a a -DehomBasis for N. 
Then {Vf''^, V^""^, 1//^^} is a a-Gröbner basis of M. 

Proof. The claim follows from the chain 

where the second equality follows from Proposition [71 (3). □ 

Lemma 12 Let N be a graded submodule of F and let 14, . . . , 14 G A^ be 
non-zero homogeneous vectors. Then the following conditions are equivalent. 

(1) The set {Vi, . . . ,Vt} is a a -DehomBasis for N . 

(2) The set {(LT^(14))'^^\ . . . , (LT^(14))^^*} generates (LT^(Ar))^^t . 

Proof. By Proposition [Hl (3) the set {(LT^(14))'''S • • • , (LT^(14))'''*} coincides 
with {LT^(14<i^^),...,LT^(14'^^'^)}, and by Proposition [1(3), (5) we have the 
equality (LT^(A^))"^* = (LT^(1/'i"'^) 1 1/ G A^, 1/ homogeneous) . The conclu- 
sion follows immediately. □ 

Theorem 13 Let N be a graded submodule of F , let Vi, . . . ,Vt & N be non- 
zero homogeneous vectors, and let us consider the following conditions. 

(1) The set {14, . • • , 14} is a a-Gröbner basis of N . 

(2) The set {14, . . . , 14} is a a-SatGBasis for N . 

(3) The set {14, . . . , 14} is a a -DehomBasis for N . 

(4) The set {14^^*, . . . , 14^^^*} generates N'^^ . 

Then we have the following chain of implications. 

(1) =^ (2) =^ (3) =^ (4) 
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Proof. The implication (1) =^ (2) is obvious. 

To prove (2) =^ (3) let v G N"^^^ . By Proposition [91 (6) and the assump- 
tion, we have v = V'^'^^ with V & N. Then there exists an index i such 
that LT^iVi) I LT^{V). Consequently LT^iV^"^) \ LT^{V'^''^), and the proof 
is complete. 

To prove (3) =^ (4) we use the equivalent condition of Lemma [12] and pro- 
ceed by contradiction. Let U G A^*^*** be a homogeneous element with minimal 
{LTa{U)Y''^^ among the elements in A^*^*** and not in (Vf*^*, . . . , V/^*). We ob- 
serve that LT^(f/) G (LT^(A^'^^*)) C (LT^(iV))'^'^* and therefore, by assump- 
tion, there exists i such that {LT-sr(^Ui)y^^ divides LT-s^{U) . We deduce that, 
for suitable c & K and t G the vector V = U — ctUi has the proper- 
ties: V G A^'^^* ; V i {Vl''\ . . . , V/^*) ; LT^(V) <^ LT^([/) . By Definition [Hl it 
foUows that (LT^(r))'^^* <^ (LT^(f/))'^^*, a contradiction. □ 

In the next example we show that the implications of Theorem [13] cannot be 
reversed. 

Example 14 Let P = Q[x,y,z], a = Lex. We use a single homogenizing 
indeterminate which we call h and we write P = Q[x,y,z,h] according to 
Remark [H then a = DegLex. Let Fi = xK^ — -^2 = x^h — , and let 
be the ideal of P generated by {^1,^2}. If F3 = y^h^ — it is easy to 
check that F3 G A^ and (LT^(Ar))«'^t = ((LT^(Fi))'^'^*, (LT^(F3))^^*) = {x.y^). 
Lemma Wl\ implies that {^1,^3} is a cr-DehomBasis for A^; however, it is not 
a a-Gröbner basis of any module, therefore (3)^^(2). Moreover, it is easy to 
see that Fi = Ff*, F2 = Ff-^ that (^1,^2) = A^*^^*, but {Ff^'^.Ff^} is not 
a (T-Gröbner basis of A^'^'^^. Therefore (4)^^(3). 

Now let P = Q[x, y,z], a = DegRevLex and let P = Q[x, y, z, h] . In this case 
we have W = DegRevLex (see Remark [1]). Let Fi = x"^ — yh, F2 = xy — zh, 
let A^ be the ideal of P generated by {^1,^2}, and let F3 = y'^h — xzh, 
so that Fl*^* = y^ — xz . Then {Fi,F2,F3} is the reduced Gröbner basis 
of A^, while {Fi, F2, Fl'^*} is the reduced Gröbner basis of a module A^ such 
that W""^ = N^^'K Therefore (2)^(1). 

We are going to use the above results to produce a strategy for computing 
Gröbner bases. First, we introducé a definition. 

Definition 15 Let a be a module ordering on T"(ei, . . . , e^) , and let W be 
its extension by W. Let G = {vi, . . . ,Vs} be a set of non-zero elements in F 
(respectively in F), and let u,v be elements in F (respectively in F). Then 
u is said to be a remainder of f by G if it is the output of the division 
algorithm applied to v and G . In that case we write u = Rem(u, G) . 
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A reordering of the elements of G may produce different elements which can 
be called Rem(ï;, G) (see for instance (llll ). Theorem 1.6.4), and a variant of 
the division algorithm, which reduces only the leading terms, may also produce 
other elements. Therefore Rem(ï;, G) would really be a set. However, for the 
sake of simphcity we write u = Rem(f , G) instead of u G Rem(f , G) to mean 
any remainder with the property that LTo-(n) not divisible by LTo-(fj) for 
all Vi e G. 



2.2 Self Saturation 



Now we write a general version of Buchberger's Algorithm. Instead of using the 



stepwise description given in the books (llll ) and (Il2l ). we prefer to concentrate 
on the main ingredients. In this way it will be easier for the reader to under- 
stand the variations presented below. We recall the notion of S- vector S(u, v) 
of M, V (see (fï?), Definition 2.5.1). If LM,, (m) = CytuGi and LMo-(t;) = c^L^Ci, 
then S(m, w) = u — ]EI}ihjÈAy_ if jj^ y g^Ye homogeneous vectors, some 



observations on the S-vector S{Ü, V) are contained in (Il2l ). Remark 4.5.3. 



Theorem 16 (Body of Buchberger's Algorithm) Let ui, . . . ,Us be non- 
zero vectors in F (homogeneous non-zero vectors in F ) and let M be the 
submodule of F (graded submodule of F ) generated by {ui, . . . , Us} . 

(1) (Initialization) Pairs = 0, the pairs; Gens = {ui, . . . ,Us) , the genera- 
tors of M ; 

Q = (/} , the a-Gröbner basis (a -Gröbner basis) of M under construction. 

(2) (Main loop) While Gens ^ and Pairs ^ ^ do 
(2a) choose w G Gens and remove it from Gens, 

or a pair {vi,Vj) G Pairs, remove it from Pairs, and let w = S{vi,Vj) ; 
(2b) compute a remainder v := Rem(w,^),- 

(2c) if V O add v to Q and the pairs {{v,Vi) \ Vi E Q} to Pairs. 

(3) (Output) Return Q. 

This is an algorithm which returns a a-Gröbner basis (a-Gröbner basis) 
of M , whatever choices are made in step (2a) and whatever remainder is 
computed in step (2b). 



Definition 17 Let ^ be a tuple of homogeneous vectors in F and V a 
homogeneous vector in F . 

(1) We call weak saturating remainder of V with respect to ^ a vector 
obtained in the following way. At each step of the division algorithm, 
the remainder is substituted by an element with the same saturation. We 
denote it by WeakSatRem(V, ^) . 
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(2) We call saturating remainder of V with respect to Q , and denote it 
by SatRem(V, , a vector (Rem(V, ö))'^^* . 



Now we describe useful variants of Buchberger's Algorithm. 

Definition 18 Let Ui, . . . ,Us be homogeneous vectors in F and let be 
the graded submodule of F generated by {Ui, . . . , Us} . If step (2b) in Buch- 
berger's Algorithm is replaced by 

(2b') compute V := WeakSatRem(Vr, ö) ; 
the procedure is called a Weak Self-Saturating Buchberger's Algorithm 
(WeakSelfSatBA). And, in particular, if it is replaced by the following spe- 
cial case of (2b') 

(2b") compute V := SatRem{W, G) ; 
the procedure is called the Self-Saturating Buchberger's Algorithm 
(SelfSatBA). 

A motivation for these names comes from the following result. 

Theorem 19 Let Ui, . . . ,Us be homogeneous vectors in F and let N be the 
graded submodule of F generated by {Ui, . . . ,Us} . 

(1) Every WeakSelfSatBA applied to {Ui, . . . , Ug) computes a a -DehomBasis 
for N . 

(2) SelfSatBA applied to {Ui, . . . , Us) computes a a-SatGBasis for N . 

Proof. To prove (1) note that, when we substitute a vector with another 
with the same saturation, the two vectors have the same dehomogenization. 
This implies that every reduction V := WeakSatRem(PF, ^) mirrors a re- 
duction of W^*^^ by Q'^^^ = {Lf^^^ \ U E Q} with only one possible excep- 
tion: though V'^^^ might still be reducible by Q'^'^^ , we may choose not to 
substitute V with an element with the same saturation (which would al- 
low the "mirror" reduction by ö), we go to step (2c) and add V to Q . In 
this case the "mirror" reduction will be later performed as a pair. Note that 
since LT(l^) is not divisible by any leading term in Q this process terminates 
by Dickson's Lemma, and the output is a set of vectors {Vi, . . . , Vt} such that 
{Vf . . . , 1//"'^} is a a-Gröbner basis of (f/f^^, . . . , Uf""^) which is A^^^"^ by 
Proposit ion [71(3). To prove (2) we observe that all the replacements of Rem 
with SatRem are equivalent to having added some element V^^^ to Gens and 
having chosen it in step (2a) just before choosing V which would consequently 
reduce to O via . □ 

Remark 20 Reconsider Fi = xh'^ — z"^ ■, F3 = V^h^ — in P = Q[x, y, z, h] , 
a = DegLex from Example [TH and run WeakSelfSatBA. 



11 



• In (2a) we choose W = Fi, in (2b') we get V = Fi, and in (2c) we add it 

to g. 

• In (2a) we choose W = F^, in (2b') we get V = F^, and in (2c) we add it 
to Q and (^1,^3) to Pairs. 

• In (2a) we choose W = S{Fi, F3) = xz^ — y^z^h, in (2b') we have these two 
reduction steps: Wi = Wh^ - Fiz^ = -y^z^h^ + z^; W2 = Wi + F^z^ = O 
and we are done. 

The output is {^1,^3} which is not a a-Gröbner basis (see Example [H]) . 
We are ready to state the main result in this section. 

Theorem 21 Let Vi, . . . ,Vs be non-zero vectors in F , let M be the submod- 
ule generated by the set {vi, . . . ,Vs} , and let {Vi, V2, . . . , VJ} be the output 
of any WeakSelfSatBA appHed to the set {v^""^ , v^""^ . . . , v^""^} . Then the 
set {Vf'^^, . . . , l^'^'''^} is a a-Gröbner basis of M . 

Proof. Let N = {v^"''', v^""" v^""^) and let {Vi, . . . , Vt} be the output of 
a WeakSelfSatBA algorithm applied to {-y]^""^, . . . , w^""^} . Theorem [JÜ] 
implies that the set {Vi,...,Vt} is a cr-DehomBasis for A^, i.e. that the 
set {Vf''^, V^^^^, . . . , is a a-Gröbner basis of A^^^^. The conclusion fol- 

lows from Proposition [TT] □ 



3 The Sugar Strategy 

If we look at the variants of Buchberger's Algorithm (see Definition [T8i) . we 
note that they differ from the ordinary algorithm (see Theorem [Tü]) only be- 
cause they allow the replacement of a vector with another one with the same 
saturation. Such replacements may create vectors with a different degree, and 
hence the corresponding critical pairs and reductions have also different de- 
gree. We observe that a reduction can also be viewed as a special S-vector 
as shown in the proof of Proposition [231 so we can concentrate on S- vectors. 
The idea is that we want to keep the original degree every time we actually 
perform such a replacement. Now it is time to become formal. 

Definition 22 Let V , V' he homogeneous vectors in F. Then V' is said to 
be a companion vector of V if there exist non-negative integers Si, . . . ,Sm 
such that V' = yl' ■ ■ ■ y^V . 

Proposition 23 Let be given a variant of Buchberger's Algorithm. For each 
homogeneous vector V which is used during the execution of the algorithm, 



12 



there exists a unique companion vector V^™ (here means sweetenedj 
which obeys the following rules. 

(1) For every input vector Ui, . . . ,Us we have Uf^ = Ui. 

(2) For every pair of vectors U, V we have S{U, Vf" = S(?7™, V^"") . 

(3) During the execution of the algorithm, when a vector V is suhstituted 
by another vector V' with the property that [V'Y^^ = V^^^ , if we have 

= yl^ ■ ■ ■y!^V'''^^ , V' = ■ ■ ■y'^V''^ with suitable non-negative in- 
tegers ai, . . . , a^, bi, . . . ,bm, then we have {V'Y^ = u'i ■ ■ ■ y'^V^^^ where 
) = Top ((ai, . . . , a^), (&i, . . . , bm)) ■ 

Proof. We need to prove that for each creation of a new vector during the exe- 
cution of the algorithm, a unique companion vector is defined. This statement 
is true for the input vectors by (1) and for the S- vectors by (2). Replacement 
of a vector with another one with the same saturation is taken care by (3). 
Every step of reduction is of the type U — ct'V with c & K, t' G T™"''". It 
can be viewed as S{U, V) and has been considered in (2). □ 

Definition 24 We denote the degree deg■^F{V^'") by sugar(V") , and we denote 



the degree deg^(S(V;, V,)^™) = deg^(^(\/r, y/™)) by sugar(\/,, \/,) . We say 



that we use the sugar strategy if the choice of the pairs in step (2a) is made 
starting with the lowest sugar, not the lowest degree. 

Elementary properties of the sugar are contained in the following proposition 
which turns out to be particularly useful for a good implementation. 

Proposition 25 Let be given a variant of Buchberger's Algorithm and let 
U,V & F be homogeneous non-zero vectors which are used during the execution 
of the algorithm. 

(1) For every U we have (Jj^^y^^ = JJ^^^ and sugar(f/) is componentwise 
greater than or equal to deg^^{U) . 



(2) Suppose that U is reducible by V , let LT-^{U) = t' Ul-^iV) , with t' = 
VT ■■■y'^r^t, t e T". Let LC^(f/) = c LC^(r) and let A = U -ct'V be 
the result of the reduction. Then we have the equality 



Proof. Property (1) foUows as an immediate consequence of Definition [22| so 
let us prove property (2). Let = y{^ ■ --yl^U , = vT ■ ■ 'Vt^V ■ Then 



SU; 



.gar(A) = Top (sugar(f/), deg^(t) + sugar(\/)) 



we have 



m 



t LT^iV) 



(1) 
(2) 
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Moreover, 



sugar(A) 

= sugar(t/ - ct'V) = sugar(5([/, V)) = deg^ {S{U, Vy^) 
= deg^ (SiW-, Vn) = deg^ ( Icm ( LT^(f/^™), LT^iVn)) 

Using formulas (1) and (2) we get 

= Icm (2/ï^+"^ ■ ■ ■ yl;;:'-^- t UT^iV),yl' 
= Icm (|/[^+"^ ■ ■ ■ 1/:^+"'" t LT^iV),yl' 

Consequently 

sugar(74) 

= deg^ ( Icm ■ ■ ■ ï/:;r+'^™ t LT^(y), 1/r ■■■y'-t LT^{V))) 

= Top ((ri + ai,...,rm + am)+ deg-^^(t V), (si, ...,Sm)+ degp^(t V)^ 

= Top (sugar([/), deg|^(t) + sugar(V)) 

where the last equality follows from formulas (1) and (2). □ 

Example 26 Consider the polynomial ring P = K[yi,y2,xi,X2] graded by 
W = (q ° J • Let U = yly2xl — y\x2 , V = y2Xi — X2 . We observe that U is 
homogeneous of degree (4,1) and V is homogeneous of degree (1,1). With 
c = 1, t' = yfxi , t = Xi we have the reduction A = U—yfxiV = yfxiX2—yfx2 
which is homogeneous of degree (4, 1). Now we consider two cases. 

Case 1 Assume that = U, 1/™ = yiV so that sugar([/) = (4,1) and 
sugar(V^) = (2, 1). According to Proposition [251 (2) . we have sugar(74) = 
Top ((4,1), (1,0) + (2,1)) = (4,1). 

Case 2 Assume instead that f/™ = U , = yiy2V so that we have sugar(?7) = 
(4,1) and sugar(l^) = (2,2). Then LT^(f/^^) = yly2xl, LT^(l^"^) = 
yiylxi. Their fundamental syzygy is {y2,—yiXi) whose degree is (4,2). 
This fact is in agreement with Proposition [251(2) for which sugar(yl) = 
Top ((4,1), (1,0) + (2,2)) = (4,2). It is interesting to observe that a 
rule of the type sugar(t V") = deg()f:) + sugar(V^) would have lead to 
sugar(?/^xiV") = (5,2), wrongly suggesting that the sugar of A should 
have to be (5, 2) . 



■ ■■i/^:rLT^(r)) 

■■■2/-tLT^(V)) 
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4 Single Gradings 



In this section we restrict our attention to the case of positive N-gradings 
i.e. gradings defined by a row matrix W witli positive entries. Then we have 
a single homogenizing indeterminate which will be called just y . A first con- 
sideration in this direction was made in Remark Hl but we can say more. 

Theorem 27 Let W G Mati,„(Z) be a row matrix with positive entries and 
let P be graded by W . 

(1) If V is a non-zero vector in F , we have LT-^{v^°'^) = UT^{v) . 

(2) If N is a graded submodule of F , and Q = {Vi, V2, . . . , Vt} is a homo- 
geneous a-Gröbner basis of N, then the set {V^'^*, V^*^*** . . . , V^'^'^*} is a 
a -Gröbner basis of N^^^ . 

Proof. Claim (1) is clear. To prove claim (2) we let be a vector in A^*^^*; 
we need to show that UTaiY^^^) \ lllaiy) for some i G {1, . . . ,t}. Proposi- 
tion[9l(4) implies that y"" - V & N for some a G N. As a consequence y"" ■ 
LT^(y) G LT^(iV) , hence there exists 1^ G ^ such that LT^(Ki) | y^'-UT^iV). 
Now, y)( LT^(l^/'^*) by (a) applied to ï; = V^^^^^^, hence LT^(\//^*) | LT^(V') 
and this concludes the proof. □ 

Corollary 28 Let N be a graded submodule of F , and let {Vi,V2, ... ,Vt} be 
the output o/ SelfSatBA applied to a set of homogeneous generators of N . 

(1) We have Vi = Vf^^ for i = 1, . . . ,t . 

(2) The set {Vi, 1^2, • • • , Vt} is a a -Gröbner basis of N""^^ . 

Proof. Using Thereom[T9]we deduce that {Vi, V2, . . . , Vj} is a a-Gröbner basis 
of a graded submodule N of F such that N^^^ = N^^^ . On the other hand, 
by construction SelfSatBA produces as output saturated vectors. Therefore 
Vi = Vf^^ ioT i = 1, . . . ,t. Now we can use the above theorem to deduce that 
{^1; V2, . . . , Vt} is a homogeneous a-Gröbner basis of N^^^ = N^'^^ , and the 
proof is complete. □ 

Corollary 29 Let vi,...,Vs be non-zero vectors in F, let M be the sub- 
module generated by {vi,...,Vs}, and let {Vi,V2, . . . ,Vt} be the output of 
SelfSatBA applied to the set {v^""^, v^""^ v^""^} . Then {Vi, V2, . . . , V} is a 
W -Gröbner basis of M'^"™. 

Proof. It follows from Corollary [28l (2) and Proposition [9l(7). □ 

Example 30 The following example shows that the above theorem and its 
corollary cannot be extended to N'" -gradings defined by matrices with m > 1. 
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The main reason is that (1) of the above theorem is not true anymore. Let 
P = K[xi,X2, X3] with 

a = Ord (^0 1^ and W = {l ï) 
If we let P — K[yi, 1/2, Xi,X2, x^] , we have 

(1 1 1 1 \ 
1 ï 1 
O O O O 1 I 
00010/ 

We considcr the ideal / generated by {xiXs — yiy2X3, y2xl — yiXi} in P. We 
check that the element x|x3 — y^xa is not in 7, but the element ^2(2^2^3" 2/1 ^3) 
which is equal to X3{y2xl — yiXi) + yi^xix^ — yiy2Xz) is in I and therefore 
the element x\x-i, — yfx^ is in P^^ . Consequently, if we let v = x\ — xi, 
we see that = y2x\ — yiXi and hence LT^(-;;'^°™) 7^ LTcr(f). Moreover 

{xiXs — yiy2X3, y2xl — yiXi} is the reduced a-Gröbner basis of / and both 
polynomials are saturated, but it cannot be the reduced CT-Gröbner basis 
of P^^ , since we have just seen that I ^ P^^ . 



5 Strategies and Timings 

In this paper we restrict our investigation and implementation in CoChA. to 
the case of the single grading. The implementation is prototypical and it is 
planned to include its final form in the forthcoming CoCbA. 5. 

We have already mentioned that a way to compute a Gröbner basis with 
inhomogeneous data is to homogenize the input data, compute the Gröbner 
basis and then dehomogenize the result. This strategy is achieved using a Weak 
Self-Saturating Buchberger's Algorithm where the choice is to never saturate 
and choose the pair or generator of lowest degree in step (2a). 

For degree compatible orderings and inhomogeneous input, the Self-Saturating 
Buchberger's Algorithm is nothing but the Standard Buchberger's Algorithm 
with sugar. In step (2a) we choose the pair or generator of lowest sugar. The 
usage of homogeneous data makes the computation of the sugar slightly more 
complicated. The result is a small overhead. 

Even if we said that we always saturate, we do not need to saturate after every 
reduction step, but we saturate only at the end when the vector (or polyno- 
mial) is no longer reducible, thus avoiding the costly operation of saturating. 
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The file containing the text of the examples discussed here can be found at 



http : //cocoa . dima . unige . it/research/papers/BigCabRob09 . cocoa 



The cl example is the classical cyclic 7 system, non homogeneous. Examples 
mora9, hairer2, Butcher and Kinl are well known in the hterature. Exam- 
ple t51 is an implicitazion problem. Example Lex is a zero dimensional ideal 
in a polynomial ring with three indeterminates, whose Gröbner basis is com- 
puted with respect to Lex, while Elim is an elimination problem with three 
polynomials in five indeterminates. 



A and H stay for the sugar and homogeneous version of the Standard Buch- 
berger algorithm respectively and S for the self-saturating version. For every 
example we examine some experimental data about the Buchberger's Algo- 
rithm performance, namely cardinality of a reduced Gröbner basis (before 
dehomogenizing in the H and S cases), the number of S-polynomials reduced 
and the number of pairs considered during a run, plus the time spent during 
the computation. The timings are for a special version of the CoCoALib-0.99 
on a Intel Core2Duo system with 2MB RAM running Linux openSUSE 10.3. 



c7 A H S hairer2 A H 



GBLen 209 443 209 72 506 72 

PolyRed 2060 2199 2060 560 3149 560 

Pairslns 61549 97910 61549 6905 127765 6905 

Time 4.52s 3.18s 4.76s 1.26s 16.00s 1.40s 



t51P 



A H 



GBLen 
PolyRed 
Pairslns 
Time 



6 80 58 
76 239 191 
242 3160 1715 
3.40s 5.23s 0.64s 



mora9 



A 



H 



2266 3977 2552 
22099 43513 25350 
2657075 7906276 3368371 
9.91s 33.05S 12.64S 
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Butch A H S 



Kinl A H S 



GBLen 
PolyRed 

Pairsins 
Time 



23 188 28 
369 987 516 

5635 17578 9256 
1.85s 3.38s 1.40s 



43 477 135 
625 4418 1471 

10779 113526 23124 
3.34s 10.80S 1.57s 



Lex 



A 



H 



Elim 



A 



H 



GBLen 
PolyRed 
Pairsins 
Time 



4 122 122 
409 345 345 
1465 7381 7381 

4.40s 0.59s 0.61s 



99 353 353 
845 1488 1488 
14330 62128 62128 

68.80S 24.45S 24.56s 



The first two Gröbncr bases are computed with respect to the DegRevLex 
ordering; we notice that the self-saturating algorithm behavior is the same 
as the Standard algorithm, the only difference being some overhead in the 
saturating case, due to more complex sugar computations, as expected. 

The last four Gröbner bases are computed with respect to lexicographic or 
elimination orderings; we see that in several cases the saturating algorithm 
offers an efficiënt alternative to the standard/homogenizing algorithm. 

Wc also notice that the numbcrs cntcrcd into the tables offer a vcry partial 
indication of the complexity of Gröbner basis computations: only in two out 
of six cases the fastest algorithm is the one with the lowest indicators. In fact, 
the possibility of using a large pool of reducers can lead to faster reductions 
and hence to an overall better performance of the algorithm itself, as it is 
clearly seen in most of the H and S verslons. 

Quantity has a quality all üs own. 
(losif Vissarionovich Stalin) 
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cations. In the "realm" of Bruno Buchberger one cannot refrain from trying 
to improve Buchberger's Algorithm. 
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